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Abstract: We carry out the renormalization and the Symanzik 0(a)-improvement pro- 
gramme for the static vector current in quenched lattice QCD. The scale independent 
ratio of the renormalization constants of the static vector and axial currents is obtained 
non-perturbatively from an axial Ward identity with Wilson-type light quarks and various 
lattice discretizations of the static action. The improvement coefficients Cy^^ and by^^ are 
obtained up to 0((7g)-terms by enforcing improvement conditions respectively on the ax- 
ial Ward identity and a three-point correlator of the static vector current. A comparison 
between the non-perturbative estimates and the corresponding one-loop results shows a 
non- negligible effect of the 0(5Q)-terms on the improvement coefficients but a good accu- 
racy of the perturbative description of the ratio of the renormalization constants. 
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1. Introduction 



Semileptonic decays of S-mesons constitute a very important source of experimental in- 
formation in i?-physics. They have been and are currently being investigated as a part of 
the research programmes of BaBar [1] and CLEO [2]. The prototype for such decays is 
—5- 7r~£"'"i/. Once the amplitude of this process is known, an experimental measurement 
of its branching ratio allows in principle to extract the CKM matrix element \Vub\- From 
a theoretical point of view, the transition is mediated by the heavy-light vector current, 
and the problem of knowing the decay amplitude amounts to calculating the QCD matrix 
element 

{n{p)\V,\B{k)) = (^k + p-q^^^^^^ /+(g^) + "^^^~"^' g,/o(g^), (1.1) 

or equivalently the form factors f+/Q{q^)-, with q = k—p the 4- momentum transferred from 
the i?-meson to the pion. 



Given the large mass of the 6-quark, a direct lattice calculation of Eq. (1.1) requires 
tiny lattice spacings (a <^ l/(5GeV)) and big volumes (L > 1.5 fm) in order to correctly 
reproduce the quark dynamics without squeezing the 5-meson at reasonably small light- 
quark masses. Various solutions have been proposed to overcome this difficulty: the reader 
is referred to [3,4] for recent reviews. Among these we mention the Heavy Quark Effective 
Theory (HQET) and the Step Scaling Method (SSM). 

In HQET Eq. (|lTD is expanded in inverse powers of the 5-quark mass. The leading 
contribution, also known as the static approximation, describes the heavy quark in terms 
of a renormalizable effective field theory. Although lattice simulations in the original for- 
mulation [5] were hampered by large statistical fluctuations due to self-energy effects of the 
heavy propagator, thanks to the recent introduction of new lattice regularizations [6] it is 
now possible to simulate static quarks with much improved numerical precision. 

The SSM, a relativistic technique based on finite size scaling, has been proposed some 
years ago by the Tor Vergata group in relation to a study of the heavy-light decay con- 
stants [7] and meson masses [8]. It has been recently shown in [9] that combining the SSM 
with HQET enables a strict control of the mass extrapolations and a consequent reduction 
of the corresponding systematic uncertainties. 

From this point of view it would be of considerable interest to extend the combined 



approach "HQET + SSM" to Eq. (ll.l|) , since a first attempt to apply the Tor Vergata 
method to the form factors has been recently presented in [10]. The goal is ambitious in 
that observables such as Eq. (|l.l| ) are intrinsically more complex than a decay constant or 
a meson mass, owing to the appearance of an additional mass scale to be identified with q'^. 

In this paper we concentrate on HQET. In view of a non-perturbative computation 
of Eq. ( |l.lD , the static vector current must first be non-perturbatively renormalized. This 
task has been partially accomplished, since in the static approximation the spatial compo- 
nents V^*'^* are renormalized by the same renormalization constant Z^*^* as the temporal 
component of the static axial current A^^^. The Renormalization Group (RG) running 
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of the latter has been computed both in the quenched approximation [11] and with two 
dynamical quarks [12]. In order to compute the renormalization constant Zy^^ of the tem- 
poral component of the static vector current Vq**^*, we derive an axial Ward identity (WI), 
much in the spirit of [13,14], relating Z^^^ to Zf^^. The scale independent ratio Z^^^/Z^^^ 
is then computed through an explicit implementation of the WI in the Schrodinger func- 
tional (SF) at the chiral point. On-shell 0(a)-improvement at zero light-quark mass is 
obtained by adding a single counter-term to the static vector current, proportional to the 
improvement coefficient c^^^, which is then tuned according to the request that the axial 
WI be satisfied at finite lattice spacing up to 0(a^)-terms. 

The improvement of the static vector current V^^*^* at non-zero light-quark mass, re- 
alized in principle through the introduction of a second improvement coefficient 6^^*, is 
not easily achievable in terms of the WI, which takes its simplest form in the chiral limit. 
For this reason, we adopt a different improvement condition, i.e we obtain by'^* by im- 
posing that the ratio of a three-point SF correlator of the static vector current at zero 
and non-zero light-quark mass be the same in two different static regularizations up to 
0(a^)-terms, thus determining the difference A6y^* corresponding to the chosen actions. 
This procedure repeats the one adopted in [6] for the determination of In order to 

isolate the value of 6^^* corresponding to all the static discretizations, the knowledge of 
^stat jg required for at least one of them. This is a difficult problem, which we solve only 
approximately by computing at one-loop order in perturbation theory for the static 
actions with the simplest lattice Feynman rules, i.e. the Eichten-Hill (EH) and the APE 
ones. This somewhat unsatisfactory solution introduces 0((7o) systematic uncertainties, 
which are discussed in detail. 

Other appealing applications where the static vector current plays a role can be found 
within the domain of twisted mass QCD [15,16], where the static axial current acquires 
a vector component after a twist rotation of the light-quark fields. This is the case, for 
instance, with the computation of B^^^, for which the matrix elements of the AB = 2 
four-fermion operators have to be normalized by appropriate bilinear correlators of the 
static axial current [17]. 

The paper is organized as follows. The axial WI is derived in sect. 2, where the no- 
tation is also established. Its implementation in the framework of the SF is discussed in 
sect. 3. Sect. 4 is devoted to a one-loop perturbative analysis of the lattice artefacts in var- 
ious WI topologies. In sect. 5 we present our non-perturbative results for the improvement 
coefficient c^*^* and the 0(a)-improved ratio Z^^^ /Z^^^, and in sect. 6 we discuss the im- 
provement coefficient b^^^. Conclusions are drawn in sect. 7. Additional tables containing 
perturbative and non-perturbative results have been collected in appendix A. 

2. Formal derivation of the axial WI 

As for the theoretical derivation of the axial WI, we follow the approach of [18,19]. For 
the moment no attention is paid to the specific regularization of the theory. We assume 
a fermion content with an isospin doublet of degenerate light-quarks ip'^ = {ipi,'ip2) and 
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a single heavy-quark, described by a pair of static fields (V'hjV'h)- order to set up 
the notation, we introduce the light-quark isovector axial and vector currents and the 
pseudoscalar density 

A^ix) = Vi(x)7M75^rV(x) , (2.1) 
V^ix) = ^Pix)J,^T'^^ix) , (2.2) 
P»(x) = Vi(x)75Vv(x) , (2.3) 

as well as their heavy-light companions (for which we explicitly indicate light-quark flavour 
indices) 

Af(x)=4(x)7/.75V'h(^), (2.4) 

V/:''{x)=Mx)7^M^), (2.5) 

P^\x) = Mxh5Mx) , k = l,2. (2.6) 

The general WI follows from the invariance of the path integral representation of 
correlation functions under chiral rotations of the light-quark fields. In particular, we 
consider an axial variation 

6l^P{x) = u;'^(x)ir"75^(x) , 6li;{x) = o;»(x)V;(a;)75^r'^ , (2.7) 

where denotes a Pauli matrix acting on the isospin space and u;"(x) is a smooth function 
which vanishes outside some bounded region R. Since the Pauli matrices are traceless, 
the functional integration measure is invariant under Eq. ( |2.7D and we conclude that the 
correlation function of a given operator O satisfies the equation 

{061S) = {510) , (2.8) 

where 

51S= [ dS cj''(x){-9^^"(x) + 2mP"(x)} (2.9) 

represents the axial variation of the light-quark action and m denotes the PC AC quark 
mass. We assume in what follows that O factorizes into the product of two operators Oi^t 
and Oext) polynomials in the basic fields and localized in the interior and exterior of R 
respectively. Accordingly, Eq. ( |2.8D reads 

(OintOext<5^5) = (Oext^^Oint) • (2.10) 

We now concentrate on the isovector component a = 1. In our specific application we 
choose Oint{x) = Vq^{x) for a; € -R and C'ext(y) = P^'^iu) for y ^ R, thus obtaining 

{Al\x)P''\y)) =2{V,^\x)P''Hy) [ d'z{d,Al-2mP'}) . (2.11) 

JR 
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If we further require -R to be a time-oriented cylinder with periodic b.c. in space, i.e. 



R = {x: ti<xo< t2}, (2.12) 
we immediately see that the space derivatives of the light axial current on the right hand 



side of Eq. (2.11) drop out, while the temporal derivative gives rise to a boundary contri- 



bution. After a space integration of both sides over x, we arrive at our final expression 

(Qf (xo)P^2(2/)) = 2{Q},\xo)P''\y) | [Q'M - Qiih)] - 2m^dSpi(z)|) , (2.13) 
where xq G [^1,^2], Ho ^ [^11^2] and we have introduced the axial and vector charges 

QUxo) = J d^x Ag(x) , (2.14) 
Qf (xo) = I d^x Atix) , (2.15) 
Qf (xo) = I d^x Fo'^x) . (2.16) 



Eq. ( 2.13| ) has to be understood as a relation among renormalized quantities. It should be 



observed that Q\ and consist of two contributions in the flavour space, corresponding 
to the non-zero matrix elements of r^. Out of them, only those with flavour content tJj2''Pi 
contribute to the right hand side of the WI. These will be denoted respectively Qj^ and P^^. 

3. Lattice implementation in the SF 

The axial WI admits a straightforward lattice implementation. We adopt here a SF topol- 
ogy where periodic boundary conditions (up to a phase 6 for the light-quark fields) are 
set up on the spatial directions and Dirichlet boundary conditions are imposed on time at 
Xo = 0, T. For a discussion of the original application of the SF to the simplest WI, namely 
the PCAC, we refer the reader to [20]. Unexplained notation closely follows [21]. 

Although the SF is formally defined in the continuum, we find it convenient to work 
at finite lattice spacing. Light quarks are assumed to be described by the 0(a)-improved 
Wilson action, with the usual Sheikholeslami-Wohlert term in the bulk and boundary 
counter-terms proportional to the improvement coefficients ct — 1 and ct — 1. No background 
field is assumed. The static quark is instead defined in terms of the action 

X 

where the forward and backward covariant derivatives 

dW^(x) = i [l^(x,0)V'(x + aO) - V(x)] , 

D^*ip{x) = i iIj{x) - W^'ix - aO, 0)V'(x - aO) , (3.2) 
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depend upon a parallel transporter W, which can be variously defined. In this paper we 
consider four possible versions, namely EH, APE, HYPl and HYP2, respectively corre- 
sponding to 

W^^{x,0) = U{x,0) , 
W^^^{x,0) =V{x,0) , 

W^^^Hx,0)=Vr''{xM,=io.7,,o.e,o.3)^ 

^"^''(-'0)=^r"(-'0)|.=(i.o,i.o,o.5)- (3-3) 

In the above definitions V{x, 0) represents the average of the six staples surrounding the 
gauge link U{x,0), while V^^^{x,0) denotes the temporal HYP link of [22], with the 
approximate SU{3) projection of [6]. 

In order to translate Eq. ( |2.13| ) into the language of the SF, we insert the static vector 
current in the middle of the bulk, i.e. at xq = T/2. The support region R is then defined by 
localizing ti and t2 at different points, with the understanding that < ti < xq < ^2 < ^ in 
order to avoid possible contact terms. The pseudoscalar density is replaced by a boundary 
source uniformly distributed over the spatial coordinates, i.e. 

^'' = IjE^h(y)75C2(z) . (3.4) 

On-shell 0(a)-improvement of the quark currents requires the introduction of operator 
counter-terms, whose structure has been discussed in [20,23]. Accordingly, we introduce 
the 0(a)-improved currents 

Af\x) = Al\x) + acj,6Al'{x) , 6Ai^{x) = '^{do + d*o)MxM,{x) ; (3.5) 

A^'^-^x) = Al\x) + acTSAl\x) , SA^ix) = M^hn^^^J + V*)Vh(x) ; (3.6) 

C^'(^) = Vo^'^ix) + 5V^o''(^) , = M^hjli^j + V*)^h(x) ; (3.7) 

and the 0(a)-improved charges 

Qf\xo) = a^^Af\x) = Q'i{xo)+acM{xo) , (3.8) 

X 

QT''\xo) =a^Yl ^o'^'(^) = Qfi^o) + ac^rSQTixo) , (3.9) 

X 

Qf ^^(xo) =a^Yl ^o''^'(^) = (^o) + a4*^*5Qf (xo) , (3.10) 

X 

where, as also explained at the end of last section, the notation O^^ always refers to a 
flavour content ipiipj- The improvement coefficients Ca, cf^^^ and c^^^ depend on the gauge 
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Figure 1: Diagrammatic representation of the SF correlation functions of Eq. ( 3.13 ). A single 
(double) line describes the propagation of a light (static) quark. 



coupling and are perturbatively expanded according to 

c = c«5o' + c('^5o' + 0(5o') ■ (3.11) 

In view of phenomenological applications, it is useful to allow for renormalized currents 
at non-zero light-quark mass. 0(a)-improvement requires in this case the introduction of 
additional mass counter-terms, proportional to rriq = m — rricr- The relations between 
renormalized currents and their bare counterparts explicitly read 

^oIrV) = -^a[1 + ht^am^A'^'^{x) , 

A'^^{x) = Zr [1 + bT'am^]A^''-'\x) , 

Fotf (x) = Z^'-'il + b'^^'am^]V,'''\x) . (3.12) 

The SF implementation of the axial WI is then realized through the introduction of a 
set of two- and three-point correlation functions, 

C(^o,yo) = 1^ E {Qv''\^o)P'Hy)^''') , (3.13) 

y 

which are graphically represented by the Feynman diagrams of Fig. |l[ It should be observed 
that the two-point correlator h\ satisfies the relation h\ = _2y^*^*'^ with f^^^'^ defined 
in Eqs. (3.22-3.24) of [23]. Once the renormalized currents are expressed in terms of the 
bare ones, the axial WI takes the form of a constraining relation among renormalization 
constants. In the chiral limit it reduces to 

h\{x,) - Z^tatZ, + ^ • ^^-'^^ 

In order to pursue a numerical implementation of Eq. (|3.14D , some geometrical pa- 
rameters have to be specified, namely the ratios T/L, ti/T, t2/T and the 0-angle of the 
SF. Concerning the latter, we consider three possible values, i.e. = 0.0,0.5,1.0. The 
other parameters will instead be collectively referred to as the topology T of the WI. In 
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Table |T] we list four possibilities. Each of them affects the noise-to-signal ratio of the non- 
perturbative simulations in its own way and introduces specific cutoff effects in the ratio 
TZ at finite lattice spacing. Therefore, a convenient choice of T imposes - at least theoret- 
ically - a balance between the minimization of the lattice artefacts and the maximization 
of the numerical signal. 



T 


T/L 


xo/T 


h/T 


t2/T 


1 


1 


1/2 


1/4 


3/4 


2 


2 


1/2 


1/4 


3/4 


3 


3/2 


1/2 


1/3 


2/3 


4 


3 


1/2 


1/3 


2/3 



Table 1: Some topologies T of the Wl. 



We remark that Eq. ( 3.14| ), which we use in order to determine Cy^*, depends as well on 



the improvement coefficients Ca and c^*^*. These have been determined respectively in [24] 
and [6] and are taken as input parameters here. In particular, c'^^^ is analytically known 
at one-loop order for the EH and APE actions, and effectively up to 0((7Q)-terms for the 
HYPl and HYP2 actions. Scaling tests of cj'^* have been extensively discussed in [6], to 
which the reader is referred for details. Here we stress that the lack of a full knowledge of 
^tat introduces systematic uncertainties at order 0(5q) in the determination of Cy*^*. On 
the other hand, the WI is independent of the boundary improvement coefficients ct and ct- 
This has been explicitely checked in perturbation theory. 

4. One-loop perturbative analysis of the WI 

A first indication of the cutoff effects related to a given choice of the topology T can be 
obtained in principle from a one-loop perturbative calculation of the WI. We anticipate 
that once the 0(a)-improvement has been carried out, the residual lattice artefacts of O(a^) 
have comparable size in the various topologies, so that a conclusive argument for the choice 
of the preferred T has to follow from non-perturbative considerations. To show this, we 
expand the ratio TZ in powers of the coupling, i.e. 

7^ = 7^w + 5o'7^(^) + o(5o') . (4.i) 

Each term of the perturbative expansion is a function of the bare quark mass m and must 
be computed at m = rricv Since the latter depends in turn upon the bare coupling, each 
correlator h of Eq. ( ^.13 ) has to be expanded according to 



h = /iWU=o + go W'^ + + /il'^l + 0{gl) , (4.2 

L J m=0 
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Figure 2: Continuum approach of 7?.''°-' with topology T — 2. 

where dm indicates a partial derivative with respect to m and the subscript "b" denotes the 
contribution of the boundary counter-terms proportional to ct — 1. The one- loop critical 
mass m-cJ^ is defined here by requesting that the 0(a)-improved PC AC quark mass vanish. 
Its values at finite lattice spacing are taken from [25,26]. 

The ratio TZ is expected to be tree-level improved, since all the improvement counter- 
terms start at 0(5(q). This expectation is confirmed by Fig. where the approach of 7^^'^^ 
to the continuum limit is reported for the topology T = 2. We observe that the slope of 

7^(o) 

increases with 6 (the scaling is perfect aX 6 = 0.0) and is independent of T. It follows 
that, in order to identify a better T, at least the one-loop contribution has to be worked 
out explicitly. 

The perturbative expansion of the two-point correlator /i^ has been discussed in [23] 
and will not be reviewed here. The one-loop coefficient of the three-point correlator /lyp 
receives contributions from Feynman diagrams corresponding to self-energy and tadpole 
corrections of the single quark legs, plus vertex corrections with gluons propagating from 
one leg to another. 

stat fl") 

Several possible improvement conditions may be imposed in order to tune Cy so 
that the 0(a)-improvement is realized at one-loop order. After some attempts, we found 
that a reasonable definition is to enforce the equation 

7^W(0l,a/L) = 7^W(02,a/i) + O[(a/L)2] , (4.3) 

which defines Cv*^*^^^ up to 0(a/L)-terms. Cutoff effects with the WI topology T = 2 are 
reported in Tables |6[]^ and in Fig. |3| for the EH (left plot) and APE (right plot) actions 
and three possible choices of the pair (0i, ^2)- 
The other topologies show similar lattice artefacts. As 
expected, different definitions converge to the same 
continuum limit, which is very small for the EH dis- 
cretization, if compared to the size of the cutoff effects, 
and somewhat larger for the APE action. It follows Table 2: Cv''*''^'' for the EH and APE 
that the extrapolation of the lattice points to the con- actions, 
tinuum is difficult and one should not expect a high 

numerical precision. In order to reduce the size of the lattice artefacts, we have employed 
the blocking procedures of [27,28]. Results are reported in Table ^. Our determination of 
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Figure 3: Continuum approach of c 



stat(l) 



for the EH (left plot) and APE (right plot) actions. 



Plots refer to the topology T = 2. Different choices of the 6'-angles provide independent definitions 



of 



stat(l) 



Plotted points correspond to L/a = 8, ... , 32. 



Cv,E*H^^ is in good agreement with the original estimate given by [29] in the framework of 
NRQCD. 

Once the improvement coefficient Cv*^*^^^ is known, the ratio TZ^^^ can be calculated 
in the 0(a)-improved theory. In Fig. ^ its continuum approach is plotted vs. (a/L)^ for 
all the topologies and the 0-angles. Corresponding data are reported in Tables |8[^ The 
main feature of the plots is the similarity of the various definitions, which differ by just a 
few percent at the coarsest lattices. Nevertheless, some topologies are more sensitive to a 
change of 9 than others, e.g. T = 4 looks almost fiat at 9 = 1.0, while it has the largest 
slope at 9 = 0.0. Remarkably, the 
spread between different T's almost van- 
ishes around 9 = 0.5, thus suggesting 
that this 0-value could be the most sta- 
ble against variations of the topology 
beyond perturbation theory. 

We extract the common continuum limit of TZ^^^ via the afore-mentioned blocking tech- 





EH 


APE 


^stat 1 ^sta± 


l- 0.0521(l)g2 


1 - 0.0093(2)gg 



Table 3: Zf ^'/Z^*'^* at one-loop order for the EH and 
APE actions. 



niques. Our best estimates are IZ 



(1) 



0.0644(1) and U 



(1) 



0.1072(2). In order to isolate 



the ratio of the static renormalization constants Z^^^^ / Z^^^ , the one-loop contribution of 
-0.116458 [30,31], has to be subtracted from 7^^^^ Results are reported in 



Za, i.e. Z, 



(1) 



Table |3[ The value obtained with the EH action is not novel: it checks the one previously 
found in [14,32] within 3%. 



5. Non-perturbative determination of Zf'^^/Z- 



stat 
V 



In order to simulate the WI non-perturbatively, we first address the choice of the geo- 
metrical parameters. Some numerical attempts suggest that the topology with the best 
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Figure 4: Continuum approach of 7?,*^-'^) for the EH (upper pfots) and APE (fower plots) actions. 
Plots refer to various SF topologies and 9 angles. Plotted points correspond to L/ a — 10, ... ,32 
for r = 2, 4 and L/a = 12, . . . , 32 for T = 1, 3. 



signal-to- noise ratio is also the one with the smallest aspect ratio T/L. This is largely 
expected on the basis of [6], since the loss of signal is mainly related to the temporal ex- 
tension of the static propagator, which for every T goes from the boundary to xq = T/2. 
Given the lack of a clear indication from perturbation theory concerning the preeminence 
of a specific topology over the others, we decide to just follow the criterion of the signal-to- 
noise ratio and to consequently adopt T = 1 for our non-perturbative study. Simulation 
parameters are collected in Table They have been taken from [6] and correspond to a 
physical size L = 2Linax = 1.436ro of the SF. Moreover, Csw is non-perturbatively tuned 
according to [24] and the boundary improvement coefficients ct and Ct are respectively set 
to their two- and one-loop values [21]. 

It is worth noting that with T = 1 the WI can be simulated directly at the chiral point, 
with no need for a mass extrapolation in the way of [33]. The k- values which have been 
used are the ones reported on the left of the third column and correspond to Kcr obtained 
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L/a [3 k 9 

8 6.0219 0.135081,0.1344011 0.0,0.5,1.0 

10 6.1628 0.135647,0.1351239 0.0,0.5,1.0 

12 6.2885 0.135750,0.1353237 0.0,0.5,1.0 

16 6.4956 0.135593,0.1352809 0.0,0.5,1.0 

Table 4: Simulation parameters used for the non-perturbative study of the improvement coefheients 

^stat^ ^stat ^j^g j.g^^iQ ^stat/^stat_ 

from the 0(a)-improved PCAC relation [11,34]. 

We also observe that the simulation at /? = 6.1628 cannot be actually performed with 
ti = r/4 and t2 = 3T/4, since these are non-integer multiples of the lattice spacing in this 
particular case. To avoid the problem, we take here ti = 2a and t2 = 7a. This choice is 
theoretically sound since no contact term turns up in the WI. It amounts to changing the 
definition of c^^^ by an 0(a)-term and of the improved ratio Zf^^'/Z^^*' by an 0(a^)-term 
at that given /?. Other choices are possible as well. The present one has the a posteriori 
advantage that it makes the /3-dependence of Z^'^^/Z^'^^ smoother than other definitions. 

To achieve a non-perturbative estimate of Z^'^^ / Zy^^ we first have to properly tune the 
improvement coefficient Cy^*. We follow the perturbative definition introduced in Eq. ( [4.3| ) 
and impose the improvement condition 

7^(^l,/3) = 7^(e2,/3) + 0(a2) , (5.1) 

where, as already explained in sect. 3, the coefficients c^*^* and Ca are taken as input 
parameters. We stress once more that since 0^^^* is known up to 0((7g)-terms, this introduces 
a systematic uncertainty in our computation, thus making the numerical estimate of c^^^ 
only non-perturbatively effective. In principle, it could be possible to avoid this by enforcing 
two simultaneous conditions, i.e. 

71(^1, /3) =71(02, /3) + 0(a2) , 

n{eup) = 7^(03, + o{a^) , Oi^e^^e^, (5.2) 

from which c^^* and cS^^^ could be determined at the same time with no approximation. 
Unfortunately, the resulting expressions for the improvement coefficients are quite involved 
and characterized by a very poor signal. For this reason we are forced to resort to Eq. (^.1|). 
Being TZ linearly dependent on c^^*, i.e. 

lZ{e,f3) = r{e,f3) + cT\l3)s{e,(3) , (5.3) 

obtain the improvement coefficient from the equation 



we 

_stat 



. - Z^'^-^^'^l + 0(a) . (5.4) 

5(6*1, /3) - s(6'2,/5) 
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Figure 5: /3-dependence of Cy'^* for some choices of the pair (6*1, 62) and of the static action. For the 
sake of readability, points (diamonds) have been slightly shifted along the horizontal axis. Dashed 
curves represent quadratic fits. 



Results at the simulation points are reported in Table of appendix A. Statistical errors 
have been computed through the jackknife method. In Fig. ^ the /3-dependence of Cy*^* is 
shown for different choices of the angles (^1, 62) and for different static discretizations. The 
most noticeable feature seems to be the large discrepancy with respect to the perturbative 
estimates given in the previous section. We also observe that the EH determination is quite 
distinct from the other regularizations, which are instead close to each other. 

The numerical values of c^^* should be independent of the choice of {61,62) up to 0(a)- 



effects. Therefore, the difference Ad 



,stat 
V 



_stat I 



-V l(6»i,6»2) 



_stat I 



is expected to decrease 



at larger /3-values. This is confirmed by our data. 



As Table 12 shows, the improvement condition with the best signal-to-noise ratio 
is the one corresponding to (^1,^2) = (0.5,1.0). This is also the one with the smallest 
perturbative cutoff effects. A quadratic fit of it in the range of the Monte Carlo simulations 
(6.0 < (3 < 6.5) leads to the parametrization 



stat 



stat 
^V,APE 



0.694 - 0.732X + 0.330a;^ , 
0.421 - 0.531X + 0.360a;2 , 



(5.5) 
(5.6) 
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Figure 6: 0(a)-improved ratio of the renormalization constants Z^'^*' / Z^^'^*' for various static 
actions. Different choices of 6 correspond to independent definitions of the WI. For the sake of 
readabihty, points (diamonds and triangles) have been sHghtly shifted along the horizontal axis. 
Dashed curves represent quadratic fits. 
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Cv,HYPi = 0.453 - 0.584X + 0.421x2 , (5.7) 
Cv,hVp2 = 0.494 - 0.528X + OAOAx'^ ; x = /3 - 6 . (5.8) 

The 0(a)-improved ratio of the renormalization constants Z^^^/Z^^^ corresponding 
to this choice of c^^* is shown in Fig. ^; the same data are collected in Table ^ To 
extract this ratio out of TZ, we used the ALPHA determination of Za reported in [35]. 
Since now all the improvement counter-terms have been taken into account, the definition 
of Zf^^/Zy^^ with 6 = 0.0 has to agree up to 0(a2)-terms with those at 6* = 0.5 and 
6 = 1.0, which have been used in order to tune the improvement coefficient. Indeed, it 



can be seen from Table 13 that the differences are zero within the statistical errors. Aside 
the non-perturbative determination, also the one-loop estimates of Table |3| are reported in 
Fig. ^ The agreement is good with EH static fermions in the whole region explored by the 
Monte Carlo simulations. It is good as well with the APE action at the largest /^-values. 
A quadratic fit (in the range 6.0 < (3 < 6.5) gives 

[Z^'^V^v'^Ieh (fo) = 0.953 + 0.0417X - 0.0828x2 , (5.9) 

[Zf''V^v''lAPE(5o) = 0.958 + 0.113X- 0.126x2 , (5.10) 



-13- 



[Zl*^V^nHYPi(5o) = 0.963 + 0.109X- 0.131x2 , (5.11) 
^^stat/^statj^^^^ (5o) = 0.961 + 0.129X - 0.146x2 ; ^ = /3 - 6 . (5.12) 



Eqs. (5.9- 5.12 ) reproduce the numbers of Table 13 within the statistical errors. For con 



venience, we also report a parametrization of /Z^ with all the improvement coeffi- 
cients set to their respective values, but c^^^^ = c^^^ = 0. With this choice, the WI is not 
O(o)-improved. Definitions corresponding to different choices of the 0-angle differ now by 
0(a)-terms, which are well above the statistical uncertainties. At 6 = 0.5 we find 

^^stat/^statj^^ (^^^^ = 0.818 + 0.186X- 0.134x2 , (5.13) 

^^stat/^stat]^^^ (5o) = 0.917 + 0.176X - 0.165x2 , (5.14) 

[Zi*^VZ^*^*]^^,pj5o) = 0.926 + 0.178X- 0.177x2 , (5.15) 

[^stat/^statj^^^^ (5o) = 0.978 + 0.158X - 0.179x2 ; X = /3 - 6 . (5.16) 



,stat 
A • 



Our final results, represented by Eqs. (p.9|)-( [5?T2 ) depend upon the choice of c 
Adopting the determination of [6] introduces a systematic uncertainty at 0((7q), which 
propagates to the ratios and can be easily estimated. We first observe that the 

improvement coefficient Cy^^ is very sensitive to variations of cj'^*. This is not surprising, 
since both the counter-terms proportional to c^*^* and cf"^* are meant to cancel the 0(a)- 
lattice artefacts of the WI. Therefore, a change of 0(1) in cf^^* is expected to produce a 
variation of the same order in d^'^^ via Eq. (|5.1|) . In practice, setting c^*^* = lowers the 



estimates of Eqs. (5^)-(5^) by 30% at the coarsest lattice spacing. Nevertheless, if the 
new values of c?, I introduced in the WI and the counter-term of the static axial 

current is explicitly dropped out at denominator of Eq. (3.14), the variation of Z^^^ /Z^'^^ 



amounts to at most 1%. This is due to a very large numerical cancellation of cf^^^ within 
the ratio TZ. It follows that Eqs. ( |5.9| )-( |5T^ ) can be assigned a systematic uncertainty of 
1%. 

Concerning the systematic uncertainty of c^^^, which has a strong dependence upon 
i^stat Q^j. determination, we naively expect that physical matrix elements of the static 
vector current will be only slightly affected by variations of Cy^^, in strict analogy with 
the static axial current, where a change of the operator counter-term is compensated by 
an opposite variation in the renormalization constant, as shown in [36]. Unfortunately, we 
have no quantitative elements at the moment to clarify if this is the case also for the static 
vector current. 



6. The improvement coefRcient b 



,stat 
V 



The axial WI at non-zero light-quark mass is complicated by the presence of a mass term 
proportional to the temporal integral of the SF correlator hyp introduced in Eq. ( 3.13| ). 



Since the integration region covers the whole interval [ti,t2]) an integrable contact term 
raises at yo = xq = T/2. Managing the integral can be disadvantageous in some cases: for 
instance, in perturbation theory it requires a complete one-loop calculation for each value 
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of the integration variable, since no Fourier transform is defined in the SF along the time 
direction. Therefore, in order to improve the static vector current out of the chiral limit, 
it is easier to look for some more comfortable observable. 



One attractive possibility is to consider a three-point SF correlator 
with the insertion of the static vector current in the bulk. To this aim 
we define 

M\x^,m) = {Y.'^W^''\x)Y}^^) , (6.1) 



) and 



70 



75 



75 



where has been introduced in Eq. (|3] 

^Ec2(y')75c;(z') 



v21 



Figure 7: Three- 
point SF correlatoi; 



jiator. 
(6.2) 



is a relativistic pseudoscalar boundary source localized at xq = T. Since we are interested 
in massive light-quarks, we keep the mass dependence explicit in the definition of M^. 
The flavour structure of the chosen valence operators allows for just one Wick contrac- 
tion, depicted in Fig. |^. To have it renormalized, all the logarithmic divergences must be 
subtracted, both those related to the static vector current and the ones induced by the 
boundary sources. In the 0(a)-improved theory the renormalized correlator reads 



M^(xo, ma) = Z^{1 + h(^am^f{l + 5^*^*amq}M^(xo, m) . (6.3) 

In order to get rid of the renormalization constants, we construct the ratio of at 
two different values of the renormalized light-quark mass, i.e. Lrn-a = 0.24 and Lm^ = 0. 
This is not sufficient to isolate 6^^^*, since the improvement of the boundary light-quark 
source contains h(^, which does not drop out in the ratio. Nevertheless, 5^ is independent 
of the static action. Therefore, it cancels when we enforce the improvement condition that 
the ratio of the three-point SF correlator be the same with two different static actions Si 
and S2 up to 0(a^)-terms, i.e. 



^ +'^^'^^"^^Mi(r/2,mer) 



. M^(T 2,m) 
1 -r v,S2 "^JmI r/2,mcr 



+ 0(a2) . (6.4) 



In the above equation, we decided to place the operator insertion in the middle of the 
bulk and to choose T/L = 1. This improvement condition provides a non-perturbative 
definition of A6^*^* = - . 

Simulations have been performed according to the parameters reported in Table ^. In 
particular, the K-values on the right of the third column correspond to Lm^ = 0.24, with 
m^ the quark mass renormalized in the SF scheme at scale // = l/(1.436ro). Numerical 
results of A^^'^* are reported in Table |l^ for the various independent combinations of the 
static actions and the usual 0-angles. They are also represented in Fig. |8|, where Ab^'^^/gQ 
is plotted vs. g^. A remarkable feature of the results is their flatness in g^. Since 

br = l + bT^'^9l + Oig'o), (6.5) 



-15- 



EH-APE I □ I APE-HYP1 i— Q- 
EH-HYP1 I O I APE-HYP2 i-^ 
EH-HYP2 I I HYP1-HYP2 i— A- 



< 



0.3 
0.2 
0.1 





I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I 



6 = 1.0 A. ^ A 



I .... I .... I .... I .... I .... I .... I .... I .... I .... I .... I 



0.3 h 

C\i o 

-5) 0.2 



0.1 




6 = 0.5 



=E©= 



0.9 0.91 0.92 0.93 0.94 0.95 0.96 0.97 0.98 0.99 1 1.01 1.02 



I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I 



■ ... I .... I .... I .... I .... I .... I .... I .... I .... I .... I .... I 



0.9 0.91 0.92 0.93 0.94 0.95 0.96 0.97 0.98 0.99 1 1.01 1.02 



0.3 
0.2 
0.1 




6 = 0.0 



-mr- 



I .... I 



— e— - 

I .... I 



I ' ' ' ' I 

=E© 



— e— -- 



I .... I 



0.9 0.91 0.92 0.93 0.94 0.95 0.96 0.97 0.98 0.99 1 1.01 1.02 

Figure 8: gg-dependence of A6^*^*/(7q corresponding to various combinations of the static actions 
and for some choices of the 0-angle. To improve the readability, points (squares and upper triangles) 
have been slightly shifted along the horizontal axis. Dashed lines represent fits to a constant. 



this could be interpreted as a good signal of scaling and could lead to the prompt conclusion 
that Aby'^^ is not dominated by 0((7Q)-terms. Nevertheless, we observe that varies from 
0.924 to 0.996 in the range of the simulations, i.e. it changes by only 8%. Such a small 
variation could be well compatible with a slight change of the differences A^^*^* even in a 
region not strictly close to the scaling one. 

A second observation is that all the differences involving the EH action have a signif- 
icant dependence on 9, with spreads varying from 30% to 60% at the various bare gauge 
couplings. This is a clear indication that large non-perturbative 0(a) lattice artefacts affect 
our definition of b^^^ based on the three-point SF correlator M^. On the contrary, the 
remaining differences, involving exclusively the ALPHA actions, are much more universal 
in 6: in these cases the spread among different definitions stays always below 0.01. 

A fit of Aby^^ /qq to a constant provides an effective non-perturbative parametrization 
of the difference of the improvement coefficients in the region of the Monte Carlo simu- 
lations. Since we have no theoretical argument to privilege one particular definition over 
the others, we decide to average the results of the fits corresponding to the three ^-values 
and to assign the averages an absolute uncertainty as large as the maximal discrepancy 
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between different ^-determinations. In this way we obtain 



(A6t*^t)EH-APE = 0.19(8)<7o' , (6.6) 

(A6t^"*)EH-HYPi = 0.18(8)5^ , (6.7) 

(A6t^"*)EH-HYP2 = 0.27(9)<7g , (6.8) 

(A6t^"*)APE.HYPi = -0.004(2)<72 , (6.9) 

(A6^*"t)APE-HYP2 = 0.078(7)ff^ , (6.10) 

(A6^*^*)hypi-hyp2 = 0.082(7)5^ . (6.11) 



In order to isolate the improvement coefficient b^^^ corresponding to each static action, 
we perform an analytical one-loop perturbative calculation of 6^*^^ and &v ape- To this aim, 
we expand the three-point SF correlator in powers of gg, i.e. 



M\xo,m) = M(°)(xo,mW) + 
+ 0(<7o') , 



+ 



(6.12) 



where the perturbative coefficients of the bare quark mass m^*^-* and m^^^ are chosen ac- 
cording to Eqs. (3.29)-(3.30) of [37]. Here mR is defined as the renormalized quark mass 
in the minimal subtraction scheme on the lattice at scale /i = 1/L. In this perturbative 
calculation we impose the improvement condition 



Ml{T/2,0.2A/L) _ 



Mi(r/2,0) 



= const. + 0(a^ 



(6.13) 



with aspect ratio T/L = 1 and 9 = 0.0,0.5, 1.0. When expanded in perturbation theory, 
this equation provides for a definition of ll^^^^^^ -\- 35^^^ up to 0(a/L)-terms. Since b^^^ — 



—0.06738(4) X Cp has been previously calculated in [37], this is sufficient to isolate b' 
Lattice data are reported in Tables 10 -111 and plotted 



C ~ 

stat(l) 
V • 



in Fig. ^. Their continuum extrapolation leads to the 
estimates quoted in Table |5[ These correspond in turn 
to an exact one-loop difference 





EH 


APE 


7 stat(l) 


0.013(1) 


-0.018(1) 



stat\(l) 



0.0324(4) 



(6.14) 



Table 5 

APE actions 



^stat(l) gjj 



which is quite a bit off the central value of Eq. (|6.6| ). Clearly, this difference may be 
attributed to the presence of non-negligible 0(5Q)-terms, which in principle could be there. 
However, the systematic uncertainty which characterizes the definition of the improvement 
coefficient with the EH action prevents us from making a more precise statement. For this 
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Figure 9: Continuum approach oi b^^ ^ ' for the EH (left plot) and APE (right plot) actions. 



Different choices of the 0-angle provide independent definitions of b 
to L/a = 8, ... ,46. 



stat(l) 



Plotted points correspond 



reason, we desist from quoting a final estimate of ?>v,e*h- Instead, we use 6v,ape 
Eqs. ( |O - |6.10| ) and quote 



tstat 



tstat 



--O.014(3)5o' + O(5o') , 



0.096(8)(7o' + 0{gt) . 



to solve 

(6.15) 
(6.16) 



!,stat\(l) 



The reader should not be surprised to see that the difference (Afey^'jEH-APE given m 
Eq. ( |6.14| ) is more precise than the single values of 6y^* reported in Table |5[ Indeed, 
the continuum estimate of Eq. ( 6.14| ) has been obtained by extrapolating the difference 



of the lattice data reported in Tables |10pi| . Part of the cutoff effects drops out in this 
difference, which makes the continuum extrapolation easier. 



6.1 A scaling test for d 



,stat 



Our non-perturbative data enable a scaling test of the three-point SF correlator M^, useful 
to assess the effectiveness of our numerical determination of Cy^*. To this aim we introduce 
the ratio 

M\T/2,m^)\e, 



C(^i,^2,mR,) 



(6.17) 



M\T/2,m^)\e, ' 

which has a well defined continuum limit, with a theoretical rate of convergence propor- 



tional to 0(a ) if the light-quark action is 0(a)-improved and is properly tuned. Fig. 10 
illustrates the approach of ^ to the continuum, corresponding to the choice of parameters 
6i = 0.5, 62 = 1.0 and Lrn-H = 0.24. The left plot refers to the non-perturbative choice of 

0. 



„stat 



provided by Eqs. (5^-5^); the right plot shows the unimproved case with c 



stat 
V 



Similar plots are obtained with different 0-angles and m^. 
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Figure 10: Scaling plots for ^(0.5, 1.0, 0.24/L). To improve the readability, some of the points (di- 
amonds and triangles) have been slightly shifted along the horizontal axis. Dashed curves represent 
independent linear fits in (a/L)^. Continuum extrapolated values are also shown. 



We observe that all the static actions give comparable results and statistical uncer- 
tainties at finite lattice spacing, save for the EH one in the improved case. If we look at the 
right plot, we note that the total variation of ^ in the simulation region is only about 5%. 
This can be attributed to a significant cancellation of the 0(a) lattice artefacts between 
numerator and denominator, which on one hand gives ^ a good scaling behaviour also in 
absence of operator improvement, but on the other makes it rather insensitive to a change 
of Cy^*. Nevertheless, we find that once Cy^* is switched on, the total variation of ^ in the 
simulation region drops to 3%, corresponding to a flatter approach to the continuum. As 
it might be expected, the strongest effect of c^^^ is at L/a = 8, where the central values of 
the lattice points are shifted by 1.5-2.9%. 

7. Conclusions 

In this paper we have studied the renormalization of the static vector current and its 0(a)- 
improvement in the quenched approximation of lattice QCD. Quark degrees of freedom are 
described by lattice Wilson-type fermions in the light sector and various discretizations of 
the static fermions, including those introduced some years ago by the ALPHA Collaboration 
(APE, HYPl, HYP2). 

Owing to the chiral symmetry of the continuum theory, the RG running of the static 
vector and axial currents coincides. Since the latter has been extensively studied in the 
literature, a complete description of the renormalization factor Z^'^^ is achieved by simply 
fixing the ratio of the two renormalized currents at a given reference scale (in our study 
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/x~£ = 2Linax = 1.436ro). To this aim we make use of an appropriate axial Ward identity 
in the framework of the Schrodinger functional. The enforcement of chiral symmetry up 
to 0(a^)-terms provides us with a lever to tune the improvement coefficient Cy***. Unfor- 
tunately, the resulting determination is not fully non-perturbative, since it relies upon a 
previous computation of c^*** which is only effective, i.e. correct up to 0(5fo)-terms. With 
regard to the numerical results, a comparison of the Monte Carlo simulations with a one- 
loop perturbative calculation shows that large higher-order contributions affect c^^^ within 
the explored region of the gauge coupling (6.0 < /3 < 6.5). On the other hand, we observe 
a good agreement between the non-perturbative determination of the 0(a)-improved ratio 
^•staty'^stat g^jj^ one- loop approximation. 

The 0(a)-improvement programme is carried out at non-zero light-quark mass via 
the introduction of a second improvement coefficient 6^^*. This is tuned on the basis 
of an independent condition involving a boundary-to-boundary three-point correlator of 
the static vector current, out of the chiral limit. The coefficient 6^^* is studied at one- 
loop order in perturbation theory for the EH and APE actions. To extend our study 
to the HYP actions, where perturbation theory is not easily handled, we adopt a mixed 
strategy: the difference A6y of the improvement coefficients between two different static 
discretizations is computed non-perturbativcly and the one-loop estimate with the APE 
action is used to isolate b^^^ in the HYPl and HYP2 cases up to 0((/Q)-terms. It has to 
be said that a direct comparison of the non-perturbative estimate of (A5y^*)EH-APE with 
its one-loop value shows that the amount of such 0(g'Q)-terms could be non-negligible and 
hard to control. Nevertheless, this problem seems to characterize the EH fermions more 
than their statistically improved versions, for which a better agreement with perturbation 
theory is expected on the basis of the experience gathered by the ALPHA Collaboration 
in previous studies of the static axial current. 

Anyway, one should always keep in mind that 6^ ^* enters the improved static vector 
current accompanied by a factor of aniq, which is rather small at light-quark masses up to 
the strange one and the commonly affordable lattice spacings. In this sense, it constitutes 
a subdominant contribution, which is not expected to have a crucial effect on the scaling 
behaviour of phenomenological matrix elements of the static vector current with external 
B^- or S^-meson states. 
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A. Additional tables 



Lja 



stat(l) I 

Cv.EH |(ei,e2)=(o.o,o.5) 



stat(l) I 

Cv,EH |(ei,e2)=(o. 0,1-0) 



stat(l) I 
Cv.EH |(6(i,6(2) = (0.5,l-0) 



4 

6 

8 

10 
12 
14 
16 
18 
20 
22 
24 
26 
28 
30 
32 



-1.21318905506 X lO^i 
-8.05799533807 X 10"^ 
-5.45943895816 x 10"^ 
-4.17792968218 X 10"^ 
-3.39132673910 X 10"^ 
-2.84872496446 x 10"^ 
-2.44775851326 x 10"^ 
-2.13776226148 X 10"^ 
-1.89020434292 x 10"^ 
-1.68756672083 X 10"^ 
-1.51841738803 X lO'^ 
-1.37494665187 X 10"^ 
-1.25162477343 X 10"^ 
-1.14441713189 X 10"^ 
-1.05031230408 X 10"^ 



-5.56860354782 x 
-3.56181870915 x 
-2.14296960443 x 
-1.58212044566 x 
-1.27133821723 x 
-1.06310641447 x 
-9.08817318086 x 
-7.87813705828 x 
-6.89431955817 x 
-6.07394710919 x 
-5.37672581715 x 
-4.77520683457 x 
-4.24986678175 x 
-3.78634326450 x 
-3.37380171813 x 



10-2 
10-2 
10-2 
10-2 

10-2 
10-2 

10-3 
10-3 
10-3 

10-3 
10-3 
10-3 
10-3 
10-3 
10-3 



-1.14115660920 x 10-2 
-7.37458670533 x 10-3 
-1.11423236142 x 10-3 
-1.04666543337 X 10-* 
4.11644560339 x IQ-^ 
7.08421808901 x IQ-^ 
1.12471960839 x lO"'^ 
1.78912840345 x IQ-'^ 
2.63691059080 x IQ-'^ 
3.59366131023 x IQ-'' 
4.60431552309 x IQ-'^ 
5.63221867541 x IQ-'^ 
6.65402844941 x IQ-'^ 
7.65528806751 x 10"^^ 
8.62744511694 x IQ-'^ 



Table 6: Three different determinations of the one- loop contribution to Cy*^* with EH static 



fermions according to the improvement condition Eq. (4.3). Numbers refer to the topology T — 2 



L/a 


stat(l) 1 

Cv.APE |(ei,e2)=(o.o,o.5) 




stat(l) 1 

Cv,APE |(ei,e2)=(o.oa.o) 


stat(l) 1 

Cv.APE |(ei,e2)=(o.5,i-0) 


4 


-1.22686673290 X 10- 


1 


-6.58409712344 x 10-2 


-2.74941292579 x 10 


-2 


6 


-7.75983902179 x 10" 


2 


-3.92522157408 x 10-2 


-1.51643269042 X 10 


-2 


8 


-4.86871345865 x 10- 


2 


-2.02408469320 x 10-2 


-2.81570404300 x 10 


-3 


10 


-3.41205392729 x 10- 


2 


-1.17290094210 x 10-2 


1.82812424278 x IQ- 


3 


12 


-2.51362502419 X 10- 


2 


-6.77983138946 x 10-3 


4.26399256380 x 10" 


3 


14 


-1.89433155853 x 10" 


2 


-3.44517803597 x 10-3 


5.84345845792 x 10" 


3 


16 


-1.43758718224 X 10- 


2 


-9.98819307355 x IQ-'^ 


6.99872560335 x IQ- 


3 


18 


-1.08514987939 X 10- 


2 


8.93569798294 x lO-'' 


7.90351127450 x IQ- 


3 


20 


-8.04166342915 x 10" 


3 


2.41124455308 x 10-3 


8.64240269108 x 10" 


3 


22 


-5.74490428783 x 10- 


3 


3.66095340524 x 10-3 


9.26291571099 x IQ- 


3 


24 


-3.83000893356 x 10- 


3 


4.71111391836 X 10-3 


9.79458495765 x IQ- 


3 


26 


-2.20751253184 X 10" 


3 


5.60799879909 x 10-3 


1.02571343800 x IQ- 


2 


28 


-8.14172229101 x 10- 


4 


6.38421491915 x 10-3 


1.06644451570 x 10" 


2 


30 


3.96087085934 x IQ-* 


7.06350798921 x 10-3 


1.10266744094 x IQ- 


2 


32 


1.45762004875 x 10-3 


7.66362155507 X 10-3 


1.13514918436 x IQ- 


2 



Table 7: Three different determinations of the one-loop contribution to c^''^* with APE static 



fermions according to the improvement condition Eq. (4.3). Numbers refer to the topology T — 2 
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L/a ni^2 {0 = 0.0) n2il {o = o.s) n^^^ {e = i.o) 



4 


3.30112791641 x 


10- 


-2 


5.07890218387 x 


10- 


-2 


5.32679355672 x 


10- 


-2 


6 


3.95982301495 x 


10- 


-2 


4.76144903008 x 


10 


-2 


4.87823906356 x 


10- 


-2 


8 


5.12452416492 x 


10- 


-2 


5.53446719297 x 


10- 


-2 


5.54812561435 x 


10- 


-2 


10 


5.62747608283 x 


10- 


-2 


5.87919096424 x 


10- 


-2 


5.88023249365 x 


10- 


-2 


12 


5.88913680345 x 


10- 


-2 


6.05967926872 x 


10- 


-2 


6.05933519321 x 


10- 


-2 


14 


6.04217811204 x 


10- 


-2 


6.16508787417 x 


10- 


-2 


6.16457789162 x 


IQ- 


-2 


16 


6.13931315301 x 


10- 


-2 


6.23177973609 x 


10- 


-2 


6.23106907219 x 


10- 


-2 


18 


6.20478428063 x 


10- 


-2 


6.27659838840 x 


10- 


-2 


6.27559138013 x 


10- 


-2 


20 


6.25099352992 x 


10- 


-2 


6.30815917333 x 


10- 


-2 


6.30682137670 x 


10- 


-2 


22 


6.28481782065 x 


10- 


-2 


6.33122587177 x 


10- 


-2 


6.32956656095 x 


10- 


-2 


24 


6.31031850685 x 


10- 


-2 


6.34860186764 x 


10- 


-2 


6.34665140076 x 


10- 


-2 


26 


6.33001903502 x 


10- 


-2 


6.36202279760 x 


10- 


-2 


6.35981895849 x 


10- 


-2 


28 


6.34555419363 x 


10- 


-2 


6.37260940915 x 


10- 


-2 


6.37019044423 x 


10- 


-2 


30 


6.35802087156 x 


10- 


-2 


6.38111147867 X 


10- 


-2 


6.37851294648 x 


10- 


-2 


32 


G.3()S177180ir, x 


10- 


-2 


G.:S8Sn459()54() x 


10- 


-2 


G.. '^8529951792 x 


10- 


-2 



Table 8: Three different determinations of the one-loop contribution to the 0(a)-improved WI 
with EH static fermions. Numbers refer to the topology T = 2. 



L/a n^lL{9 = 0.0) n^lL{o = 0.5) n^lL{o = i.o) 



4 


7.83613664415 x 


10- 


-2 


9.63395380897 x 


10- 


■2 


1.02312037738 x 


10- 


■ 1 


6 


8.29999758507 x 


10- 


-2 


9.07196239420 x 


10- 


-2 


9.31211717407 x 


10- 


2 


8 


9.43090648772 x 


10- 


2 


9.79649261355 x 


10- 


2 


9.83100792082 x 


10- 


-1 


10 


9.92568717719 x 


10- 


-2 


1.01312590309 X 


10- 


-1 


1.01130674952 x 


10- 


1 


12 


1.01847718258 x 


10- 


-1 


1.03111765757 X 


10- 


-1 


1.02755357444 x 


10- 


■ 1 


14 


1.03368358402 x 


10- 


-1 


1.04185677884 x 


10- 


-1 


1.03765015807 x 


10- 


■ 1 


16 


1.04335470132 x 


10- 


-1 


1.04878533404 X 


10- 


-1 


1.04436312823 x 


10- 


-1 


18 


1.04988139034 X 


10- 


-1 


1.05352674754 x 


10- 


-1 


1.04907826792 x 


10- 


-1 


20 


1.05449158322 x 


10- 


-1 


1.05692363172 x 


10- 


-1 


1.05253903968 x 


10- 


■ 1 


22 


1.05786795092 x 


10- 


-1 


1.05944779845 x 


10- 


-1 


1.05517080769 x 


10- 


■ 1 


24 


1.06041438462 x 


10- 


-1 


1.06138003224 x 


10- 


-1 


1.05723087837 X 


10" 


■ 1 


26 


1.06238214652 x 


10- 


-1 


1.06289597523 x 


10- 


1 


1.05888244607 x 


10- 


■ 1 


28 


1.06393415756 x 


10- 


-1 


1.06411014965 x 


10- 


-1 


1.06023326229 x 


10- 


-1 


30 


1.06517980204 x 


10- 


-1 


1.06509988458 x 


10- 


-1 


1.06135695966 x 


10- 


-1 


32 


1.06619471630 x 


10- 


-1 


1.06591897793 x 


10- 


-1 


1.06230536156 x 


10- 


■ 1 



Table 9: Three different determinations of the one- loop contribution to the 0(a)-improved WI 
with APE static fermions. Numbers refer to the topology T = 2. 
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L/a 


, stat(l) / o ^\ 
Ov.EH (t' = 0.0) 






, stat(l) //) 


0.5) 






, stat(l) ,n 


1.0) 




6 


-1.92494990552 x 10 


-1 


7 


59451390152 


x 


10- 


-2 


7 


94433253811 


X 


10- 


-2 


8 


-1.28284962759 x 10 


-1 


7 


19668221779 


x 


10- 


-2 


9 


82580984223 


X 


10- 


-2 


10 


-7.83097066265 x 10 


-2 


6 


93670764950 


x 


10- 


-2 


8 


92282752595 


X 


10- 


-2 


12 


-5.20449399705 x 10 


-2 


6 


60882882658 


x 


10- 


-2 


8 


20217785512 


X 


10- 


-2 


14 


-3.59939921230 x 10 


-2 


6 


27858498787 


x 


10- 


-2 


7 


59502707254 


X 


10- 


-2 


16 


-2.52419762755 x 10 


-2 


5 


97829596381 


x 


10- 


-2 


7 


09345102910 


X 


10- 


-2 


18 


-1.76143672046 x 10 


-2 


5 


71059870738 


X 


10- 


-2 


6 


67548992369 


X 


10- 


-2 




1 1 Q8S^nfi7zl48 V 10 
1 . lyooouu ( A lu 


-2 





471 '^n81 71 '^1 


X 


10- 


-2 


L) 




X 


10- 


-2 


22 


-7.71840193796 x 10 


-3 


5 


25663880980 


X 


10- 


-2 


6 


01574097982 


X 


10- 


^2 


24 


-4.40467430491 x 10 


-3 


5 


06225291197 


X 


10- 


-2 


5 


74858429769 


X 


10- 


-2 


26 


-1.78765598057 X 10 


-3 


4 


88528309736 


X 


10- 


-2 


5 


51202560392 


X 


10- 


-2 


28 


3.08255934246 x 10~ 


4 


4 


72321849465 


X 


10- 


-2 


5 


30036608571 


X 


10- 


-2 


30 


2.00560511473 x 10" 


3 


4 


57399574759 


X 


10- 


-2 


5 


10928394946 


X 


10- 


-2 


32 


3.39251992852 x 10" 


3 


4 


43591299385 


X 


10- 


-2 


4 


93543727896 


X 


10- 


-2 


34 


4.53369168047 x 10^ 


3 


4 


30755529827 


X 


10- 


-2 


4 


77620635919 


X 


10- 


-2 



Table 10: Three different determinations of the one-loop contribution to 6^*'^* with EH static 
fermions according to the improvement condition Eq. ( 6.13| ). 



L/a 




Ca*^e^(^ = o.o) 






,stat(l) _ 

"V.APE — 


0.5) 






,stat(l) _ 

"V,APE — 


1.0) 




6 


-1 


92340408686 


X 


10- 


1 


7 


38125998350 


X 


10 


-2 


7 


32869781736 


X 


10- 


-2 


8 


-1 


36845655003 


X 


10- 


1 


6 


16212128505 


X 


10 


-2 


8 


47024088104 


X 


10- 


-2 


10 


-9 


27641005323 


X 


10- 


2 


5 


36113971906 


X 


10- 


-2 


7 


11683168091 


X 


10- 


-2 


12 


-7 


03192609445 


X 


10- 


2 


4 


68591075301 


X 


10- 


-2 


6 


11505675084 


X 


10- 


-2 


14 


-5 


67994367873 


X 


10- 


2 


4 


12586767144 


X 


10 


-2 


5 


32267457567 


X 


10- 


-2 


16 


-4 


77974262206 


X 


10- 


2 


3 


66636780399 


X 


10 


-2 


4 


69184320708 


X 


10- 


-2 


18 


-4 


14377364146 


X 


10- 


2 


3 


28289740025 


X 


10- 


-2 


4 


17877409334 


X 


10- 


-2 


20 


-3 


67710201232 


X 


10- 


2 


2 


95586046439 


X 


10 


-2 


3 


75125941270 


X 


10- 


-2 


22 


-3 


32534218174 


X 


10- 


2 


2 


67170066138 


X 


10 


-2 


3 


38739675354 


X 


10- 


-2 


24 


-3 


05477831878 


X 


10- 


2 


2 


42109755463 


X 


10 


-2 


3 


07226742530 


X 


10- 


-2 


26 


-2 


84338924925 


X 


10- 


2 


2 


19744102751 


X 


10- 


-2 


2 


79543183872 


X 


10- 


-2 


28 


-2 


67622184958 


X 


10- 


2 


1 


99587562440 


X 


10- 


-2 


2 


54935704579 


X 


10- 


-2 


30 


-2 


54281229985 


X 


10- 


2 


1 


81272438403 


X 


10 


-2 


2 


32844980335 


X 


10- 


-2 


32 


-2 


43562070602 


X 


10- 


2 


1 


64513560157 


X 


10- 


-2 


2 


12845506743 


X 


10- 


-2 


34 


-2 


34911031398 


X 


10- 


2 


1 


49085018459 


X 


10 


-2 


1 


94607725049 


X 


10- 


-2 



Table 11: Three different determinations of the one-loop contribution to 6^^* with APE static 
fermions according to the improvement condition Eq. ( 6.13| ). 
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.7 




,.stat 

X.AJ'K 


,.stal, 

WII'l'I'l 


,.stai 


(0.0,0.5) 


6.0219 

6.1628 
6.2885 
6.4956 


0.756(22) 

0.577(18) 
0.548(17) 
0.399(18) 


0.478(18) 

0.337(14) 
0.334(14) 
0.240(14) 


0.513(18) 

0.360(14) 
0.359(14) 
0.261(14) 


0.553(18) 

0.409(14) 
0.416(14) 
0.324(14) 


(0.0,1.0) 


6.0219 
6.1628 
6.2885 
6.4956 


0.707(10) 
0.566(8) 
0.532(8) 
0.406(8) 


0.433(8) 
0.328(6) 
0.318(6) 
0.242(7) 


0.467(8) 
0.352(7) 
0.342(6) 
0.263(7) 


0.508(8) 
0.402(7) 
0.398(7) 
0.327(7) 


(0.5,1.0) 


6.0219 

6.1628 
6.2885 
6.4956 


0.690(8) 
0.562(6) 
0.527(6) 
0.409(6) 


0.419(6) 

0.325(5) 
0.312(5) 
0.242(5) 


0.452(7) 

0.349(5) 
0.336(5) 
0.264(5) 


0.493(6) 

0.399(5) 
0.392(5) 
0.328(5) 



Table 12: Non-perturbative determinations of c^*'** for various gauge couplings and static actions. 
Different choices of (^1,^2) correspond to independent definitions of the improvement coefficient. 



9 


P 


r vstat / vstatl 

l^A /^v J EH 


L^A f-^v Jape 


r -zstat / -zstatl 

l^A /^V JhYPI 


r 7stat / ystatl 

l^A / JhYP2 


0.0 


6.0219 


0.9549(13) 


0.9611(13) 


0.9662(13) 


0.9643(13) 




6.1628 


0.9564(10) 


0.9725(10) 


0.9771(9) 


0.9785(10) 




6.2885 


0.9585(8) 


0.9799(9) 


0.9837(9) 


0.9859(10) 




6.4956 


0.9527(7) 


0.9823(7) 


0.9847(7) 


0.9890(6) 


0.5 


6.0219 


0.9549(13) 


0.9601(12) 


0.9651(12) 


0.9633(12) 




6.1628 


0.9562(10) 


0.9723(10) 


0.9769(10) 


0.9784(9) 




6.2885 


0.9582(8) 


0.9796(8) 


0.9834(7) 


0.9856(9) 




6.4956 


0.9528(6) 


0.9823(7) 


0.9847(6) 


0.9890(7) 


1.0 


6.0219 


0.9540(11) 


0.9601(11) 


0.9651(11) 


0.9633(11) 




6.1628 


0.9561(9) 


0.9723(10) 


0.9769(10) 


0.9784(8) 




6.2885 


0.9583(8) 


0.9796(7) 


0.9834(7) 


0.9856(9) 




6.4956 


0.9528(6) 


0.9823(6) 


0.9847(5) 


0.9890(6) 



Table 13: Non-perturbative determinations of the 0(a)-improved ratio Zf^^/Zy'^^ for various 
gauge couplings and static actions. Different choices of 6 correspond to independent definitions of 
the WI. 
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^ /3 (A6y**) EH-APE (A6y^*) BH-HYPl (A6y^*) EH-H^ 

0.0 6.0219 0.2177(14) 0.2098(15) 0.2967(18) 

6.1628 0.2122(17) 0.2091(18) 0.2925(21) 

6.2885 0.2130(22) 0.2123(22) 0.2953(26) 

6.4956 0.1869(37) 0.1890(39) 0.2626(11) 



6 /3 (A6y ^*)apE-HYP1 (A6y **) APE-HYP2 ( A6y )hyP1-H 

0.0 6.0219 -0.0079(7) 0.0787(10) 0.0865(6) 

6.1628 -0.0031(7) 0.0800(10) 0.0831(6) 

6.2885 -0.0008(8) 0.0820(12) 0.0828(7) 

6.4956 0.0021(11) 0.0756(17) 0.0735(10) 



6 (3 (Aby^*) EH-APE (A6y^*) EH-HYPl (A6y^*) EH-H-i 

0.5 6.0219 0.2094(18) 0.1987(19) 0.2826(23) 

6.1628 0.2020(24) 0.1983(25) 0.2793(28) 

6.2885 0.2050(28) 0.2018(30) 0.2828(34) 

6.4956 0.1789(51) 0.1777(52) 0.2501(57) 



^ (A6y *'')ape-hypi (A6y ^'')apE-HYP2 (A6y *'')hypi-h 

0.5 6.0219 -0.0106(9) 0.0729(13) 0.0835(7) 

6.1628 -0.0036(10) 0.0771(14) 0.0807(8) 

6.2885 -0.0032(11) 0.0776(15) 0.0808(8) 

6.4956 -0.0004(15) 0.0720(21) 0.0724(12) 



6 P (Aby"**) EH-APE (Afty"**) EH-HYPl (Afty"**) EH-H-i 

1.0 6.0219 0.1332(25) 0.1250(27) 0.2019(34) 

6.1628 0.1301(30) 0.1270(30) 0.2018(34) 

6.2885 0.1391(34) 0.1379(36) 0.2153(41) 

6.4956 0.1284(63) 0.1285(64) 0.2013(68) 



G (A6y **)ape-hypi (A6y **)ape-hyp2 (A6y **)hypi-h 

1.0 6.0219 -0.0082(12) 0.0686(19) 0.0768(10) 

6.1628 -0.0030(12) 0.0716(16) 0.0747(8) 

6.2885 -0.0012(13) 0.0761(18) 0.0772(10) 

6.4956 0.0001(17) 0.0728(24) 0.0727(13) 



Table 14: Non-perturbative determinations of Ab^^^ for various gauge couplings and static actions. 
Different choices of 6 correspond to independent improvement conditions. 
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